Shor's factorization algorithm with a single control qubit and imperfections 
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We formulate and numerically simulate the single control qubit Shor algorithm for the case of 
static imperfections induced by residual couplings between qubits. This allows us to study the 
accuracy of Shor's algorithm with respect to these imperfections using numerical simulations of 
realistic quantum computations with up to = 18 computational qubits allowing to factor numbers 
up to A*' = 205193. We confirm that the algorithm remains operational up to a critical coupling 
strength tc which drops only polynomially with logj A*'. The obtained numerical dependence of Ec 
on log2 A'^ is in a good agreement with the analytical estimates that allows to obtain the scaling for 
functionality of Shor's algorithm on realistic quantum computers with a large number of qubits. 
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I. INTRODUCTION 

Shor's factorization algorithm [i| demonstrates expo- 
nential efficiency gain compared to any known classical 
algorithm and is definitely the most important quantum 
algorithm in quantum computation Q . The possibilities 
of experimental investigations of the algorithm are rather 
restricted due to small number of experimentally avail- 
able qubits and moderate accuracy of available quantum 
gates. Thus the maximal number factorized experimen- 
tally is iV = 15 with a 7-qubit NMR-based quantum 
computer 0]. 

In view of these experimental restrictions the numeri- 
cal simulations of Shor's algorithm in presence of realistic 
imperfections becomes essentially the only tool for deter- 
mination of the conditions of algorithm operability with 
few tens of realistic qubits. The first steps have been 
done in [1, Q for factorization of N = 15. More re- 
cently, larger values of N have been studied with N up 
to 33 in and up to 247 in Q. An interesting ap- 
proach was used in [8|: the Quantum Fourier Transform 
(QFT) part of Shor's algorithm has been performed in a 
scmiclassical wa y u sing the one qubit control trick (see 
e.g. 0, [13, [ill. Il2l. [ij]) while the modular multiplica- 
tion has been performed with up to 20 qubits including 
the workspace using the circuit described in These 
works analyzed the effects of dynamical phase errors [3 
and discrete qubit flip errors [1, [l^l- Another impor- 
tant type of errors is related to static imperfections in- 
duced by residual coupling between qubits which under 
certain conditions can lead to quantum chaos melting of 
a quantum computer (isj . Such type of static imperfec- 
tions generally give a more rapid decay of the fidelity 
of quantum computation compared to random uncorre- 
lated phase errors in quantuni gates (see [l3| and Refs. 
therein). In our recent work [Tjl we studied the effects 
static imperfections for Shor's algorithm factorizing num- 
bers up iV = 943 using up to L = 30 qubits. 

In this work we combine the two approaches used in 
d, [ll] and [l3l using Shor's algorithm with a single con- 
trol qubit. This algorithm was introduced and analyzed 



in [3, [13, [m, [H, [3 ■ It allows us to perform extensive 
numerical studies of the effects of static imperfections for 
Shor's algorithm factorizing numbers up to iV = 205193 
that is significantly larger compared to [1, [3] and [l3| . 
Thus, while in [l7[ we were able to consider Uq ~ 10 
computational qubits (plus ni = 2nq = 20 control qubits 
with the total number of qubits L = 30), we use in the 
present work values up to nq = 18 computational qubits. 
Together with the single control qubit this requires a sim- 
ulation of a quantum algorithm with the total number of 
qubits L = 19. We remind that without the one con- 
trol qubit simplification this would require L = 54 qubits 
that corresponds (if simulated on a classical computer) 
to an array of 2^^ complex elements of 2^^ bytes (= 2^^ 
GB) in total (counting 16 bytes per complex double pre- 
cision number). Therefore the one control qubit simplifi- 
cation is crucial for the numerical simulation and allows 
to increase the value of iV by a factor 200 compared to 
[l7t and 10^ compared to [g|. This allows us to deter- 
mine the parametric dependence of the accuracy on the 
imperfection strength, number of qubits and number of 
gates. For this we use a simplified but generic model 
of imperfections which can be applied to various imple- 
mentations of Shor's algorithm discussed in the literature 
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The paper is structured as follows: in Section [II] we 
remind the standard Shor algorithm (with or without 
imperfections) and we explain how it is possible to ob- 
tain a modified version with a single control qubit. The 
results of the numerical simulations are presented in Sec- 
tion mil and the discussion and conclusion are given in 
Section IIVI In appendix [^ we describe the numerical 
extrapolation scheme, used in section Hill allowing to de- 
termine in an efficient way the inverse participation ratio 
for an unknown random discrete distribution. 



II. SHOR'S ALOGRITHM WITH ONE 
CONTROL QUBIT 

Let A'^ be a large integer number of which we want 
to determine its prime factors and x > 2 a small inte- 
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ger number relatively prime to N. The aim of Shor's 
algorithm [l| is to determine the period r defined as the 
minimal positive integer r such x"^ ~ 1 (modulo N). As 
pointed out in [l| the knowledge of the period allows 
(with a certain probability) to obtain a non-trivial fac- 
tor of N. The determination of r can be efficiently done 
by a quantum algorithm with Uq computational qubits, 
chosen such that N < 2"'' , and ni = 2nq control qubits. 
Actually the effective number of Uq may be larger in or- 
der to take eventual workspace qubits into account which 
could be needed to realize explicitly the quantum modu- 
lar multiplication in the computational register in terms 
of elementary one- and two-qubit gates. The actual num- 
ber of quantum gates scales like ^ (log2 A^)^, while for 
all known classical computation algorithms the scaling is 
almost exponential. However, here wc do not enter into 
these details and as in Ref. we simply assume that 
we can perform this modular multiplication operator in 
some global way not to be specified in the numerical sim- 
ulation of the quantum computation. This assumes that 
there are no quantum errors in the register with n; con- 
trol qubits. 

In order to keep the following notations simple we as- 
sociate to general operator products an ordering from left 
to right: 



n-l 



llO.li') =OoOi ...0„-i\^) 
j=o 



(1) 



This convention is necessary to keep the following de- 
scription unique and mathematically precise. 



A. Standard Shor algorithm with ni control qubits 

First, we remind the standard Shor algorithm with im- 
perfections as modeled in We start with the initial 
state 



l^o) = |0)„, |1)„, 



and then compute 

Tlj-l 



iV'i) = n i^o) 

J=0 



(2) 



(3) 
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3=0 




= t^QFT \i'2) 




ni-l ( 


C^QFT 


= ^ n 




3=0 \ 



(5) 



(6) 



Here we apply the ordering convention ([T]), Hj denotes 
the Hadamard gate acting on the j-th control qubit, 
B)jf: (ip) is the controlled two-qubit phase shift gate and 
for later use we also note the simple one-qubit phase shift 



gate as Bj (ip). The operator R reverses the order of 
the ni control qubits (see ref. p^ . section 3, for more 
notation details). Here, C/qft is the standard quantum 
Fourier transform 



The operator U. 



(3) 

Cmult V 



in Eq. Q is the controlled 



modular multiplication operator acting on the computa- 
tional register as 



\{yx) mod Af) , y==0,...A^-l 
|y> , y = iV,...,2"^-l 

(7) 

if the j-th control qubit is |1) and U^l^^^^^{x)\y) — \y) if 
the j-th control qubit is |0). The operator e"^^J denotes 
the error operator which only acts on the computational 
register (see Ref. [l7[ and the next section for details). 



The case of the standard pure Shor algorithm is simply 
obtained by putting 6Ti.j = that eliminates all errors. 

The final step of Shor's algorithm is a measurement of 
all control qubits, thus destroying lips), and resulting in 
measured numbers from each control qubit: aj G {0, 1} 
which provide the (measured) control space coordinate 
by a = X^Jlo^ "-"^i 2-' - In the pure case the probability 
distribution of a is given by 



sin^(A/fc7rar/Q) 
^ sin^(7rar/(3) 



(8) 



where Q = 2"' and Mk = [{Q-k- l)/r] + 1 ^ Q/r > 1 
only depends weakly on k. The function P{a) is com- 
posed of r well localized peaks at mQ/r with Q = 2"' 
and m = 0, 1, ...,r— 1. Since the measured value of a 
is very close to one of these peaks one obtains by a con- 
tinuous fraction expansion the value of r^rovided that 
m and r are relatively prime (see Rcfs. [l|, Il7j for more 
details). This algorithm only works with a certain prob- 
ability since m and r may have a common non-trivial 
factor or because in some rare cases even the knowledge 
of the period r is not sufficient to obtain a non-trivial 
factor of A^ In the case of imperfections {SHj ^ 0) 
the peaks of the probability distribution of the control 
space coordinate become larger and further delocalized 
secondary peaks appear. These effects of imperfections 
reduce furthermore the success probability of Shor's al- 
gorithm and can be characterized by the inverse partici- 
pation ratio which is the key quantity investigated in [T7| 
and in section IIIII of this work. 



B. Reduction to a single control qubit 

The algorithm described in Eqs. ([3]l6|) requires a large 
number ni of control qubits and it is therefore quite dif- 
ficult to implement, both in numerical simulations (on a 
classical computer) or eventually in future experimental 
realizations of quantum computers. For numerical sim- 
ulations the large qubit number is especially costly in 
terms of memory and computation time. 
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However, for this particular algorithm it is possible to 
reduce the number of control qubits to one single qubit 
using a scheme based on a semiclassical irnplementation 
of the QFT pioneered by Griffiths et al. 3 and which 
was later applied to the pure Shor algorithm by Mosca 
et al. and Parker et al. 

In order to understand this significant simplification 
we note that in the above Shor algorithm the operator 
factors associated to a particular value j in the products 
commute with the operators on the right side associated 
to j > j- This allows to regroup the operator products 
in Shor's algorithm as follows: 

n? — 1 

with operators Vj defined by: 

V, = H, J n B^^{ir2^-^)\^ (10) 

Furthermore, in Eq. ([9]) the j-th control qubit is not mod- 
ified by the later factors Vj with j < j and we can there- 
fore measure it immediately after the application of the 
factor Vj (before application of the remaining factors). 
However, after measuring this j-th control qubit, wc need 
to replace in the remaining factors Vj (with j < j) the 

(2) 

two-qubit control phase shift gates Bk^' {n2^~^) by sim- 
ple one-qubit phase shift gates which are classically con- 
trolled: B^^^\aj Tr2^~^) where aj € {0, 1} is the measure- 
ment result of the j-th control qubit. In this way, we see 
that the information obtained from measuring the con- 
trol qubit j is fed-back for use of the later values j < j 
and thus the full algorithm can be done with a single 
control qubit (with j = 0). 

Thus, from now on, we assume that the control register 
contains only a single control qubit associated to j = 0. 
This new algorithm can be put in the following recursive 
form with states \(fj) and numbers aj € {0, 1} to be 
determined as: 



|0)i|l) 







(11) 



where the first factor in \ipo) refers to the single control 
qubit. Furthermore, for j = ni — 1, ni — 2, . . . , 1, (in 
this order) we compute: 

{n,-l 
n B^,'\ak 7^2^-^^) \ X (12) 
k=j+i 



Cmult 



X 



XB^^\aj+i Tt)Hq \ipn:-j-l) ■ 

The state |(/?„j_j) is obtained from \ipni-j) by measuring 
the single control qubit and the measured value will be 



denoted by aj . Due to the projection of the measurement 
IVrii-j) bas the form: 

\^n,-j) =-\aj)\ipn,-j) (13) 

where |'^n,-j)n, is a state which only lives in the compu- 
tational register. The state will be used as the initial 
state in the next step with j — I and since aj may be 1 
the application of the Hadamard gate Hq may provide a 
"wrong sign" in this case: 

i?o|l)|^„,-,> = ^(|0)-|l)) \^n,-,) (14) 

and therefore we have introduced in Eq. (|12|1 the addi- 
tional gate B^Q^\aj tt) (for j — 1) such that: 



Bl,'\a,n)Ho\a,) 



(15) 

which is indeed the desired initial condition for the next 
step. In the above algorithm wc also introduce artifi- 
cially aj — for J = ni which is normally not relevant 
and simply provides a proper functioning of the iteration 
at the first step at j = ri; — 1. The quantum circuits 
associated to the iteration (fT2|l and the one control qubit 
Shor algorithm are shown in Figs. [T] and [5] respectively. 

The final operator i?, which inverses the order of the 
control bits, can be done classically by the reconstruction 
of the measured control space coordinate: 



"i— 1 



*ni-l- 



(16) 



The one control qubit version of Shor's algorithm repro- 
duces exactly the same probability distribution of the 
control space coordinate as the standard Shor algorithm 
described above in section Hi Al 



III. NUMERICAL RESULTS 

The probability distribution P{a) in Eq. ([8]) only 
weakly depends on k. This expression, viewed as a func- 
tion of a real variable a, has r equidistant strongly local- 
ized peaks of width unity, of height 1 jr and located at 
mQ/r with m = 0, 1, . . . , r — 1. For integer values of a 
the real peak height is probably smaller than 1/r since 
the exact position is not reached. However, the choice 
ni = 2nq ensures that there is at most only one integer 
value of a close to the exact peak [l[ . 

Wc model static imperfections generated by residual 
couplings between qubits in the frame of the generic 
quantum computer model analyzed in [isj . These resid- 
ual static imperfections produce additional unitary rota- 



•y rota- 



tions Us = e in the quantum gates. As in Refs. 
we model the static imperfections with the effective per- 
turbation operator 



5n,= Y.5.{j)a\'-^+2Y^.ma\^\. 



i+l 



(17) 
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FIG. 1: The effective gate /^(aj+i, . . . , a„,_i) which corresponds to one iteration step of Eq. (|12|l . This gate includes one 
measurement producing a classical value aj G {0, 1} and depends on the previously obtained values aj+i, . . . , an,-i. For the 
initial iteration at j = nj — 1 we furthermore put = 0. We assume that the quantum measurement of the single control 
qubit produces a normalized state obtained from a projection and the subsequent normalization. Therefore the effective gate 
conserves the normalization but is not linear and produces the classical output Oj. 
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FIG. 2: The sequence of effective gates from Fig. [T]put to- 
gether providing the one control qubit Shor algorithm and 
producing the classical binary values aj which allow to re- 
construct the measured control space coordinate according to 

Eq. (nnj. 



where ct^'''' are the Pauli operators acting on the «th qubit 
(of the computational register) and Ji{j) are ran- 

dom coefficients distributed according to: 

S^iJ), Mj) e [-V3e,V3e] . (18) 

As it was done in [l7j . we consider two models where 
the random coefficients Si{j), Ji{j) are different for each 
value of j (generic imperfection model) or equal for all 
values of j (correlated imperfection model). 

Since, for SHj ^ 0, the success of the algorithm de- 
pends essentially on the probability of hitting one of r 
peaks in the process of the a-measurement the most di- 
rect way to study this probability is by clashing all the 
peaks into one. or in other words, adding them all to- 
gether by taking a modulo s where s is the nearest integer 
value of the ratio Q/r and thus reducing all probabilities 
inside one cell with s states. In this way we obtain a new 
distribution of global search probability W{a): 

r-l 

W{a) = J2P{[a + s+jQ/r] mods) (19) 

3=0 

where now a = — s/2, . . . , s/2 — 1 (the difference of a for P 
and W is clear from the context) and s ^ Q/r is the dis- 
tance between peaks. For the ideal algorithm this global 



probability W{a) has one peak at a = while in the 
case of imperfections the main peak may become larger 
and secondary peaks appear. We also use the original 
notation of [l | putting a = c. 

As in Ref. l3| we study the "dclocalization" effects of 
quantum chaos due to the imperfections by computing 
the inverse participation ratio associated to the global 
probability distribution W{a): 



(20) 



In Ref. [13], the complete state jV^a) was calculated 
from a classical simulation (with up to 30 qubits: ni = 20 
and riq = 10) thus allowing to determine exactly all the 
key quantities such as the full probability distributions 
P{a) and W{a), inverse participation ratio, variance and 
this without actually measuring and destroying the state 
iV'a). However, the computation of these quantities is 
only possible due to the (quite expensive) classical sim- 
ulation of a guantum algorithm. In fact, the original 
Shor algorithm actually contains a measurement of the 
control space variable a. This point is rather important 
since the one control qubit version of Shor's algorithm is 
essentially based on measurements and gives a significant 
reduction of numerical computational efforts. In such a 
case, as with a real quantum computer, we are not able 
to determine directly the exact probabilities W{a) (un- 
less STij = where the theoretical formula ([8]) applies) 
but we may draw as many values a with this probability 
distribution as we want simply by repeating the simula- 
tion of Shor's single control qubit algorithm with other 
sequences of measurement outputs. Thus, many repe- 
titions of many random results of measurements is the 
prize to pay for the reduction of number of control qubits 
from 71; to one. 

Of course we may replace the exact probabilities W{a) 
by histogram probabilities which will be as accurate as we 
want provided that the number of series of measurements 
is sufficiently large. If we want to determine the full dis- 



5 



tribution the classical simulation of Shor's single control 
qubit algorithm is no longer advantageous (in computa- 
tion time) as compared to the direct simulation of the full 
Shor algorithm as done in [l7j. But if we need to know 
only certain averaged characteristics of the distribution 
W{a), e.g. the inverse participation ratio, then Shor's 
single control qubit algorithm becomes much more effi- 
cient compared to the approach used in [17[ ■ To compare 



the validity of these two approaches we verified for some 
small numbers of n; and Uq that the histogram distribu- 
tion obtained from the simulation of Shor's single control 
qubit algorithm reproduces very accurately all details of 
the exact distribution obtained from a full Shor algorithm 
simulation if both cases are simulated with the identical 
disorder realization for SHj. In Fig. [3] we show as an 
example a comparison of the distribution W(a) obtained 
from both types of simulations and for an average over 10 
disorder realizations (identical disorder realizations arc 
used for two computational methods). 




FIG. 3: (Color online) The global probability distribution 
W{c), c = a, as defined in Eq. (|19|l . averaged over Nr — 10 re- 
alizations of random static imperfections, for N = 493, x — 2, 
r = 56, e = 0.04, disorder realizations are identical for top and 
bottom panels. The top panel corresponds to the complete 
computation of Shor's algorithm as in Ht!] using in total 27 
qubits, the bottom panel corresponds to the histogram com- 
puted using A^meas ~ 12 000 measurements (error tolerance 
~ 2% for the associated inverse participation ratio) in the 
simulation of Shor's single qubit algorithm with in total 10 
qubits. The grey/green dashed line corresponds to the prob- 
ability at e = 0. 

The single qubit Shor algorithm is advantageous for the 
computation of the inverse participation ratio ^ provided 
^ is not too large because this requires less measurement 
series samples than for the full histogram to achieve a 
reasonable accuracy. However, a simple replacement of 
W{a) by the histogram probabilities in (|20p is not very 
optimal for modest sample numbers Nji since the aver- 



age of (with ^ obtained from histogram probabili- 
ties) scales with the sample number and is not identical 
with the exact value of (with ^ obtained from the 
exact probabilities W{a)). In appendix [K\ we describe a 
numerical extrapolation scheme that allows to take this 
into account and to determine a more accurate value of 
^ for a finite sample number and also to control its sta- 
tistical variance. The numerical results of ^ presented in 
the following have been obtained by this extrapolation 
scheme and choosing a sample number to achieve a 2% 
precision. 

In Fig. [3]we show the dependence of ^ on the strength 
of imperfections e. The values of ^ are obtained by the 
full simulation of Shor's algorithm as in [l3| and by the 
simulation of the single qubit Shor algorithm using the 
extrapolation scheme. We see that both methods give 
very close results even for large values of ^. Thus we 
may use the more efficient single qubit algorithm to test 
effects of imperfections for factorization of numbers N 
much larger than those of pTf . 



2000 




FIG. 4: (Color online) Averaged IPR ^ given by ([20j as a 
function of e for A*' — 493, x — 2. The full circles symbol 
correspond to the simulation of the full Shor algorithm used 
in 17] with the complete control register, therefore ritot ~ 27. 
The open circles correspond to the IPR obtained by the sin- 
gle qubit Shor algorithm with the total number of qubits 
L = 7iq + 1 — 10 and the extrapolation scheme with an error 
tolerance of 2% that required about Nmcas ~ 12000 measure- 
ments. The inset shows the dependence for small e values, 
the horizontal dashed line marks the quantum chaos border 
defined by the condition ^(e^) = 10^(0). There are Nr = 10 
disorder realizations which are identical for full and open cir- 
cles. 



The results of Figs. I3I4I show that with the increase 
of imperfections strength e the peak in the distribution 
W{a) is washed out and the quantum chaos destroys 
the operability of the algorithm. As in the criti- 
cal strength of the imperfections ec at the quantum chaos 
border can be approximately determined by the condition 
^(ec) = 10C(e = 0). In Fig. [Slwe show the dependence of 
Ec on the number N factorized by Shor's algorithm in a 
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log-log scale. The dependence on N can be described as 

_ B 

- [log, NY ^^^^ 

with numerical constants B and (3. The fit of numeri- 
cal data done for large TV values in the interval 5.5 < 
log2 iV < 18 gives B = 0.7877 ± 0.073, (3 = 1.275 ± 0.045 
for the generic imperfection model and B = 0.958±0.154, 
/? = 1.546 ± 0.08 for the correlated imperfection model. 
The values for the exponent (3 differ slightly from those 
obtained in [l3| where we had (3 = 1.420 ± 0.054 for the 
generic imperfection model and (3 = 1.523 ±0.068 for the 
correlated imperfection model. In view of strong fluc- 
tuations related to the arithmetic properties of x, r and 
N we can consider that the agreement with the results 
obtained in [l3| for not very large values of < 1000 is 
rather good. The new data allowed to increase the values 
of < 205193 by a significant factor 200 that gives more 
accurate values of the exponent (3. The obtained values 
of (3 are close to the values given by the theoretical esti- 
3| with (3=1 for the generic imperfection model 
1.5 for the correlated imperfection model. We 



mates 
and /3 

attribute the deviations of numerical values of (3 from 
the theoretical values to strong arithmetical fluctuations 
which require a large scale of log2 A^- variation. We also 
note that the statistical fluctuations related to random- 
ness and disorder in realizations of imperfections are rel- 
atively small since the standard deviation from disorder 
average gives an error bar which is approximately of the 
syinbol size in Fig. [5] (the same is true for the data of 

To give more information we present the results and 
parameters of our numerical simulations in the TableUfor 
large N values 10^ < N <2 - 10^ which were inaccessible 
in Ref. [13 • 

IV. CONCLUSION 

The extensive numerical simulations performed in this 
work allowed to analyze the accuracy and operability 
bounds for Shor's algorithm in presence of realistic static 
imperfections. The results show that above the quantum 
chaos border Cc given by Eq. [21] the algorithm becomes 
not operational while below the border the factorization 
can be performed. This border drops only polynomially 
with the logarithm of factorized number N . The alge- 
braic power (3 of this decay is close to the theoretical 
estimates obtained in [l3|. The numerical values of (3 
are close to the values obtained in [l^l where the factor- 
ization was studied for significantly smaller values of N 
compared to the present work. Due to that we think that 
our results give the real asymptotic value of the algebraic 
exponent (3 for the quantum chaos border in Shor's al- 
gorithm in presence of static imperfections. Even if the 
values of /3 = 1 or 1.5 are relatively low still the accuracy 
requirements for quantum gates become rather restric- 
tive if one wants to factorize such large N values as those 



0.2 








0.15 








0.1 










■ 






0.05 




■ Qg B i 8 




0.025 




B 


. ■ 1 





4 




6 

l<>g2 N 


8 10 


0.2 
0.15 


O 
□ 








0.1 










^ 0.05 
0.025 




□ 




□ □□ 


0.01 












4 


6 


8 10 


12 14 16 18 



log2N 



FIG. 5: Dependence of the critical imperfections strength ec 
on logj A'' in log- log scale, is obtained from the criterion 
^(ec) = 10^(e = 0), where ^ is the inverse partici pat ion ratio 
(|20p . The top panel shows the results obtained in marked 
by full symbols (full, red, circles are for the generic imperfec- 
tion model; full, green, squares are for the correlated imper- 
fection model) and the results of this work marked by open 
symbols (open circles: the generic imperfection model; open 
squares: the correlated imperfection model). The bottom 
panel shows the data of this work up to A' = 205193 (open 
symbels). The straight lines show the fits Ec = -B/(log2 A')'' 
of open symbols data for the two models of imperfections in 
the interval 5.5 < logj A^ < 18 with B = 0.7877 ± 0.073, 
/3 = 1.275 ± 0.045 for the generic imperfection model (solid, 
red, line), B = 0.958 ± 0.154, f3 = 1.546 ± 0.08 for the corre- 
lated imperfection model (dashed, green, line). 



used in classical computers (see more detailed discussion 
in [13). 

This work was supported in part by the EC IST-FET 
project EuroSQIP. For numerical simulations we used the 
codes of Quantware Library [23 |. 



APPENDIX A: EFFICIENT NUMERICAL 
DETERMINATION OF THE IPR FOR A 
DISCRETE RANDOM VARIABLE 



Let us consider a discrete random variable x with pos- 
sible values X ^ Q, 1, . . . , Q—1 and probabilities p(x) > 
properly normalized '^.j.p{x) = 1. The inverse partici- 
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TABLE I: Table for the values of ec plotted in Fig. [5] for values 
of A'^ greater than the ones of In Fig. [S] we see that below 
N « 1000 the values of quite well coincide with the ones 
computed previously in Tr|. For simplicity of the table we 
rounded the numbers to three significant digits. Number of 
realizations Nr is approximate. The number of measurements 
Nmaas givcs the Order of magnitude of the maximum number 
of measurements for both imperfections models (generic 
and correlated ^^-^ imperfection models). Here L = + 1 is 
the total number of qubits. 
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pation ratio of x is defined as: 



7 = E^'(-)- 



(Al) 



Here ^ denotes roughly the number of possible x values 
with a significant ( "maximal" ) probability. In the case 
of quantum states |?/') = i^ix) \x) withp(x) = |V'(a;)P 
the quantity ^ is also refereed as the (inverse participa- 
tion ratio) localization length. Obviously, ^ is easily com- 
puted provided the exact probabilities p{x) are known. 
In Ref. [l3|, this was indeed the case since we were able 
to calculate the full quantum state after application of 
Shor's algorithm but before measuring the control space 
variable. However, in this work, where we use the one- 
control-qubit version of this algorithm as described in 



section ini this is no longer possible since the exact val- 
ues of p{x) are not known and we are "only" able to 
draw an arbitrary number of values Xj, j = 1, . . . , i? us- 
ing this probability distribution by simply repeating the 
one-control-qubit Shor algorithm R times. In the limit 
R ^ oo this should in principle allow to recover p{x) and 
^ with sufficient accuracy but for "moderate" values of 
i? (e. g.: i? w 10^ — 100^) this is not very precise and 
can be improved by a kind of extrapolation scheme in R 
which we will now explain. 

Suppose that xi, . . . , xn are independent random vari- 
ables with the same (unknown) probability distribution 
p(x). For a given set of xi, . . . , xu (representing "nu- 
merically obtained values") we introduce the histogram 
probabilities by: 



Pr[x) 



1 

R 



1 



R 



(A2) 



where n^^^y is the number of Xj values being equal 
to X. Obviously the average of S^xj with respect to 
Xj is: {Sa:xj) = J^xjPi^j) ^xxj = p{x) and therefore 
{pr{x)) = p{x). For this simple quantity the average 
histogram value indeed coincides with the exact value. 
However, this is not the case for other quantities. Let 
us for example consider the IPR value obtained by the 
histogram probabilities: 



1 



Pr{x) 



(A3) 



with the following average: 

X 



1 




A ^ 1 



(A4) 



Here the first term arises from the (j = I)- and the second 
term from the [j ^ /)-contributions. Wee see that for 
a finite ratio i?/^ the average histogram-IPR does not 
coincide with the exact IPR. Eq. (|A4[) allows for the 
numerical extrapolation: 



Co 



l-l/R 1 
^R "i 7 — 7T, — Sfi 



1 - ^r/R 



(A5) 



where is the numerical histogram-IPR of which we 
hope that it is close to its average (for large enough R) 
thus justifying (jASp . This extrapolated IPR will be more 
reliable than the histogram-IPR for moderate values of 
R and allow for a more accurate determination of the 
functional dependence of the IPR on the different pa- 
rameters. However, is still subject to statistical er- 
rors and therefore we also need to compute the variance 
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of the histogram-IPR and related to this we also need the 
average of the second order histogram-IPR: 



42, ii ^ 



as compared to the exact second order IPR: 



(A6) 



(A7) 



This quantity is comparable to 1/^^. Actually, using 
Cauchy- Schwartz inequality (for two vectors Vx = p(a;)^/^ 
and Wx = p{x)^^^) we find: 

l^{v\v){w\w)>\{v\w)\' = ^ (A8) 

with the standard scalar product: {v\w) ~ w.^,. In 

Eq. (jA8|) we have equality if p{x) =const.> for certain 
values of x and p{x) — for the other values of x. 

Repeating the calculation (jA4|l for the second order 
IPR we find: 



with: 

h{x) = (pU^)) 

= p^pix) + 7p^il- p)p\x) + 

+6pil-p){l-2p)p\x)+ (All) 
+ {l-p){l-2p){l-3p)pHx) . 

We note that the prefactor "7" in the second term of 
(|A11[) arises from 4 permutations of the type j = I = k ^ 
m and 3 permutations of the type j = I ^ k = m \n the 
summation index. The prefactor "6" in the third term 
arises from 6 permutations of the type j ^ I ^ k ^ m. 

Furthermore for x ^ y wc obtain: 



{PR{x)PR{y)) 

p\l- p)p{x)p{y)+ (A12) 
+p(l - p)(l - 2p) [p\x)p{y) +p{x)p^{y)] 
+ (l-p)(l-2p)(l-3p)p2(x)p2(y) 



and therefore: 



1 



1 / ^ 

= E (PRix)) = ;^ E { E '^--^ '^^ 



X \j,l,k=l 



^^(i?p(x)+3i?(i?-l)p2(a 



i?3 

X 

+R{R-l){R-2)p^x] 
= p^ + 3p(l-p)- + {l-p){l-2p)- (A9) 

? ? 2 

where the three contributions correspond to the cases 
where the three values j, / and k are equal, or only two 
of them or none of them are equal. 

The evaluation of the variance of the histogram-IPR is 
more complicated but straight forward: 



(AlO) 



+4p(l-p)(l-2p) 



(A13) 



1 

e 



In the numerical scheme we determine 4/?/ and ^2 r for 
one realization oi xi, . . . , xr and using (|A4p . (|A9p we 



determine approximate values of ^ and ^2 a,nd by 
(jAlSl) the variance of ^i?/- The number R of -values 
is increased until the relative error is below a certain 
threshold, typically 2%. We note that according to (|A13p 
for the special case ^2 = the variance scales with p^ = 
R~^ and not with the usual behavior p = R~^ . Even 
for ^2 < 4^ the numerical prefactor of the R^^ term may 
be quite suppressed as compared to the i?~'^-term and 
therefore it is better to be careful and not to neglect this 
term in (|Al3jl . 
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